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This paper deals with the study of the effects of external magnetic field and gravity inclination on the chemotactic 
bacterial bioconvection by considering a Continuum model. Chemotaxis causes cells to swim out of the plume because the 
high concentration of the cells constituting the plumes leads to a lower concentration of oxygen in the surrounding fluid. 
There are cases where the layer is no longer vertical. In such cases, the gravity inclination plays a significant part. A 
similarity solution is found for the plume in which the cell flux and the volume flux could be matched to those in the 
boundary layer and also outside the suspension regions. Axisymmetric plumes are formed by applying two scales one with 
respect to the radial co-ordinate and the other with respect to the similarity variable. The effects of magnetic field and 
gravity inclination are remarkable, encouraging and the computed results are in excellent agreement with those of 
hydrodynamic case in the limiting case. 
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World's major portion consists of biomass and bioconvection is the spontaneous formation of patterns in 
suspensions of swimming microorganisms due to their tactic nature viz. oxytactic, gyrotactic, chemotactic etc., The 
microorganisms exhibiting bioconvection have the following key features: They are denser than water and they swim 
upwards due to their tactic nature. This leads to an unstable situation in the system and thus an overturning instability 
develops leading to pattern formation [1][2][3]. The study of hydromagnetic convection is of great practical interest since 
Magnetic field has a strong influence on the system in many real time situations. . 

Some literatures pertaining to bioconvection in deep chambers are [4] [5]. The study of such a phenomenon has a 
variety of applications in biological and physiological problems. Further, chemotaxis and oxygen consumption are 
important in setting up the basic state and soon after, the resulting plumes are entirely buoyancy driven and the cells are 
merely advected. In such cases, the velocity would vary across the plume [6] [7] . Experiments on bioconvection containing 
suspensions of bacteria (Bacillus Subtilis) have revealed the formation of Falling plumes when the system becomes 
unstable [8]. The present work investigates the nonlinear Hydromagnetic bioconvection to study the effect of magnetic 
field on the formation of falling plumes (Axisymmetric) where the oxygen consumption and chemotaxis are important in a 
gravity inclined environment. The model constituted the quasi - steady situation in which an upper boundary layer 
containing a high concentration of bacteria feeds a falling plume of cell-rich fluid. The suspension is divided into three 
separate regions, a cell-rich upper boundary layer of known thickness A , a falling plume of unknown width s which also 
contained a high concentration of bacteria and the fluid outside the plume which had to circulate in order to conserve mass. 
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Here, the assumption of the axisymmetric nature of the plume reduced the 3D-problem to 2D-problem [9]. Not much 
literature is available in this direction. The solutions were obtained by a Fast Computational Technique. 

2. MATHEMATICAL FORMULATION 

The bacterial suspension (Bacillus Subtilis) contained in a deep chamber reveal the development of a thin upper 
boundary layer of cell-rich saturated fluids which becomes unstable, leading to the formation of falling plumes which is a 
complex phenomenon. This was used as a basis for our mathematical model. The whole suspension is under the influence 
of uniform magnetic field. 

The dimensionless governing equations are: 

The equation of cell conservation 

SN 



dt 



= v.[h(g)vn - UN - H(e)yNve] 



(1) 



The Navier - Stokes equation (with Boussinesq approximation) 
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Equation (2) in the component form is 
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The equation of oxygen concentration 

ae 

St 



=v.[8ve-ue-H(e)6pN] 



(3) 



The conservation of mass 
V.U = 

The magnetic induction equation 

dt 



■ + 



(U.V)H = (H.V)U + P m V 2 H 



(4) 



(5) 



The variables are non - dimensionalized by using the following scales: 
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V = "^ V ' N = ^T' ° = r -r mm ' ° N = D no h ( 6 )> V = bV s H(e)V0, K = K H(9)bV s , 

iA,ui,H=i 

h »» H 

Where h : width of the chamber, V s : it has dimensions of velocity, N : initial cell concentration, U : saturated 

fluid velocity, K : the inclinatio n vector = (sin 8* ,0, cos 8* ), 5 * : gravity inclination to the vertical D N0 , K : constants, 
D N : the cell diffusivity, 9 : the oxygen concentration, C : the initial concentration, H(o) : the step function, T : the time, 
H : the constant magnetic field. 

2.1 Dimensionless Parameters 

p. W 2 ,=^,5=^,B- = ^.r= uN °»>- p ->,Sc = ^.P„ 



D c( C 0- C mi„)' D N0 D N0 Pw VD N0 ' vD N()Pw D N()' m D N0 

Where |3 : strength of oxygen consumption relative to its diffusion, y : measures the relative strength of 
directional and random swimming, 8: ratio of oxygen diffusivity to cell diffusivity, F : Bio - Rayleigh number, Sc: 
Schmidt number, P m : magnetic prandtl number, v: kinematic viscosity of the fluid, p c ,p w : densities of cell and 

water, g: acceleration due to gravity, U : volume of the cell, B * : modified Hartmann number, fx * : magnetic permeability, 
V m : magnetic viscosity of the fluid. 

2.2 Boundary Conditions 

No slip condition at Z = 1 ( bottom of the chamber) 
Stress free condition at the upper surface of the chamber i.e., at Z= 
The vertical components of velocity vanish at both the boundaries 
Zero cell - flux at both the boundaries 

Zero oxygen flux at the bottom surface and C = C at the free surface 
H = at both the boundaries 

The vertical components of velocity vanish at both the boundaries Mathematically, 
AtZ = 0, U.Z = 0, ^— (u.z)=0, 9 = 1, H(6)— -yNH(9)— = 0, (h.z)=0 (6a) 

59 dN 

AtZ=l, U.Z = 0, UxZ = 0, — = 0, — = 9, H.Z = 9 (6b) 

dZ dZ 

3. AXISYMMETRIC PLUMES USING RADIAL CO-ORDINATES 

In the plume, the radial co-ordinate is scaled as, 

R = sr, 9<S<1 (7) 
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Rescaling: N = N A n, 9 = 1 + C A C, W = W A w, U = eW A u, P = P A p 
Where N Ai C A , W A _ and P A are scale factors. 

Using (7) and (8) the axisymmetric governing equations (neglecting o(e 2 ) terms) are as follows, 
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Here, W A =8 (to retain the advection term ), C A = — (to retain the chemotaxis term ), N A = — — ( to retain 

y e 

the oxygen consumption term in 12 ), r = o(^'8 _2 ) (to retain the buoyancy term in 11), P A =e~ 4 
(to retain the pressure term in 1 1), H A = s~ 2 ( to retain the induction term in 15). 



This leads to — = , hence p = p (Z) in ( 10). 

fir 



(16) 



Substituting for C A , W A and N A in (9) and (12) we get the following equations, 
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Differentiating (11) w.r.t. r and substituting for N A W A _ H A and T we get, 
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Differentiating (15) w.r.t r and substituting for H A ,W A we get, 



Su Sw 5 2 w Sw Sw S 2 w 

1- u — — H 1- w 
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Now, imposing the boundary conditions: 

— = 0, — = 0, u = 0, = 

fir fir fir 

5C 

Also, r— > co, n— >0, >0, w->0 

fir 

3.1 Similarity Solution for Axisymmetric Case 

In order to obtain a similarity solution [ 10] [1 1] for (17), (18), (19), (20) the solution is posed in the form 

h::Z a ,w::Z b ,n::Z c ,C::Z d ,u::Z a+b+1 

(h: width of the plume, a = ¥2, b = 0, c = -1, d = 0) 



(19) 



(20) 



(21a) 



(21b) 



(22) 



Since h: Z , the similarity variable is defined as r| = — — 

Z /2 



Assuming the solution in the form 



n = Z 'H(ri), C = G(ri), i|/ = ZP(r|), u = Z 



F _ F' 



, w 



-F' 



(23) 



(24) 



( y/ : Stream functions, Primes denote differentiation w.r.t tj ) 

Substituting these into (17) and integrating once w.r.t 7 with the boundary conditions at 77 = 0, we get the 
following equation: 



HF+2r|HG'-r|H' = 0. 

Substituting into (18): 

riG" + G'-^G'F-riH = 0. 



(25) 



(26) 



Substituting into (19): 
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-F"'-^-F" + ^-F' + Sc- 
r\ r| 2 if 



Substituting into (20): 
M = F 

The boundary conditions are , 
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At r|->oo : H->0, G'->0, >0, >0 

r\ r\ 

4. SOLUTION 

Fory^O, CFD technique is employed. However for y = (i.e., when the chemotaxis is unimportant 

in the plume) analytical solutions are possible withP = 0(l), C Q =1, N D =e~ 2 and r = s 2 r. Following [4] [9] the 
solutions for the equations (23, 24, 25, 26) are found to be (see table 1) 
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Table 1: Solutions for F, H, G' 
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Also solutions satisfy the boundary conditions at r| = and r\ — »co . 



Qf oos8*(2-B*)(8-B*) 



2304 

Since j[HF'dn = -Q] 



for Sc = 1 and A 2 



QfcosS*(3-2B*)(l5-2B*) 



9216 



for Sc = 2 



In table 1, the first two columns correspond to the hydromagnetic case while the third column corresponds to the 
hydrodynamic case in the absence of gravity inclination [8] . Our results are in excellent agreement with those of the 
hydrodynamic case in the absence of gravity inclination. 
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5. RESULTS AND DISCUSSIONS 

In this study, the deep chamber experiment [8] has been modeled in three separate regions: 

a) an upper boundary layer of depth X R b) a falling plume of width e c) The region outside the plume. 

In sections 3 and 4, solutions for the cell and the oxygen concentration and the fluid velocity in the upper 
boundary layer are determined under the influence of a uniform vertical magnetic field and gravity inclination. The 

solutions are found to depend on the parameters like, Sc (Schmidt number), Q (the cell flux), T (Bio-Rayleigh number), 

B * (Magnetic Parameter), 5* (gravity inclination) and 8 (diffusivity ratio). The computations are performed using the 
MATLAB tool; the computed results are presented through graphs in Figures 1 to 13. 
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Figure 1: Fvsr/ 



Figure 2: Hvsrj 



The following observations are made: The effect of variation in the magnetic parameter B * on the profiles of 
velocity w = ^~ ^/ryj ' ^ e concen tration H and the oxygen concentration G is studied for the values, 



Sc = 5 = r = Q=1.0. Here, the oxygen concentration is considered as 9 : 



1 



where y is assumed to be 



always positive so that, all the bacteria are active. There is no effect of gravity inclination on F. In figures 2 and 3 the effect 
of gravity inclination with variation in magnetic parameter B* is studied forSc = 5 = T — Q = 1.0 and 5* =50° 

In Figure. 1, the effect of similarity variable r\ on the F profile are shown. F decreases enormously and remains 

constant as r\ — >oo for all values of B* Further, F is negative and its value is highest in the hydrodynamic case. In other 

words, F increases in absolute value as B * increases, and the vertical fluid velocity w, at the center of the plume increases 
indicating that the horizontal fluid flow into the plume increases. 

In Figure. 2 the effect of similarity variable on H profile is shown, it reveals that, as B decreases the cell 
concentration in the plume increases as expected. Physically it means that, the higher the concentration of the cells, the 
greater is the consumption of oxygen which means that the oxygen concentration at the center of the plume decreases. In 

Figure. 3 the effect of similarity variable on G' profile is shown. Clearly the width of the plume decreases as B* increases. 

The oxygen concentration is more in the hydrodynamic case ( B* = 0) when compared to the hydro magnetic case ( B* ^ 0). 
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H-* il->- 

Figure5 .Hvsrj Figure6: G'vsrj 



Figures 4, 5, 6. represent the effect of variation of T on the profiles F, H, G' for fixed values of the parameters 
Sc = 1, Q = 2, B* = 1, 5* = 50° . The values of Y considered are 0.2, 2 and 5. It is observed that there is an enhancement 

in the values of F, H, G' when compared to the uninclined case [8]. The effect of buoyancy becomes important when F is 

large. The cell concentration is more for small values of T and the plume becomes narrower for large T . When the cell 
concentration in the centre of the plume increases, the plume becomes narrower, accordingly the oxygen profile becomes 
narrower and the oxygen concentration at the center of the plume increases. The consumption of oxygen will be more. 
Therefore, the velocity of the fluid in the Center of the plume will be larger when the buoyancy force is dominant. 

But w— >0 more rapidly than for small values of T . The decrease in M for the increase in T indicates that less fluid is 
entrained by the plume. 




Figure7 :Fvs77 Figure8 : H VS77 Figure9: G'vsrj 
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Figures 7, 8, 9. present the graph of the profiles F, H, G' when the values of Sc (= 1, 2) are varied. The other 
parameters have fixed values viz., 5 = 1, Q = 1, B* = 0.2, 5* = 10° . It is observed that as in the hydrodynamic 

case (b* =o) , the variation in Sc has a significant effect on the behavior of the profiles. There is a drastic change in the 
values of F for Sc = 1 and 2. As Sc increased, F decreases rapidly and F — > a constant value as r\ — »ooas expected. The 
cell concentration will be more for Sc = 2 and accordingly the oxygen consumption in the plume will be more and there 
will be a reduction in the oxygen concentration in the plume for any value of 5* . 

Finally figures 10, 1 1 present the effect of 5* on the profiles of H and G' . Both H and G' increase with 5* & the 

general behavior is the same (as in 8*= 0). The same behavior is exhibited even in the hydrodynamic case 
(Figures 12, 13). 

Finally it is concluded that the governing dimensionless parameters viz., 5*,B*,Sc,r,5 and Q have a strong 

influence on the bioconvective system considered. It is observed that (i) the governing dimensionless parameters have a 
remarkable effect in the hydrodynamic as well as in the hydromagnetic cases (ii) the qualitative nature of the profiles is 
almost the same in both the cases but there is a drastic difference in the quantitative nature of the profiles. 




Figurel2:H vsrj Figurel3: G'vsrj 

Figures clearly indicate the strong influence of the magnetic parameter on the present bioconvective system, these 
clearly suggest that the plume convection could be suppressed or enhanced through the proper choice of the magnetic 
parameter. The results are in excellent agreement with the hydrodynamic case in the presence and absence of inclination. 
From the present study it is possible to retrieve the results pertaining to nonlinear chemotactic bioconvection for the 
following cases viz., non-linear hydrodynamic chemotactic bioconvection in the presence and absence of gravity 
inclination and non-linear hydromagnetic chemotactic bioconvection in the presence and absence of gravity inclination. 
Therefore, it is possible to suppress or enhance chemotactic bioconvection by a suitable choice of the governing 
dimensionless parameters. 
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